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Background: Liver transplantation (LT) is the only treatment option for patients with advanced liver disease. Currently, liver donation to 
these patients, considering priorities, is based on the Model for End-Stage Liver Disease (MELD). MELD score is a tool for predicting the risk 
of mortality in patients with advanced liver disease. However, few studies have so far been conducted in Iran on the efficacy of MELD score 
of these patients. 

Objectives: This study reviews the present status of the MELD score and introduces a new model for optimal prediction of the risk of 
mortality in Iranian patients with advanced liver disease. 

Patients and Methods: Data required were collected from 305 patients with advanced liver disease who enrolled in a waiting list (WL) in 
Imam Khomeini Hospital from May 2008 to May 2009. All of the patients were followed up for at least 3 years until they died or underwent 
LT. Cox regression analysis was applied to select the factors affecting their mortality. Survival curves were plotted. Wilcoxson test and 
receiver operating characteristics curves for survival predictive model were used to compare the scores. All calculations were performed 
with the SPSS (version 13.0 ) and R softwares. 

Results: During the study, 71 (233%) patients died due to liver cirrhosis and 43 (14.1%) underwent LT. Viral Hepatitis (43.7%) is the most 
common cause of end-stage liver disease among Iranian patients. A new model (NMELD) was proposed with the use of the natural 
logarithms of two blood serum variables (total bilirubin and albumin) and the patients' age (year) by applying the Cox model: 
NMELD = 10 x (0.736 x In (bilirubin) -1312 x In (albumin) + 0.025 x age + 1.776) 

Conclusions: The results of the Wilcoxon test showed that there is a significant difference between the usual MELD and our proposed 
NMELD scores (P < 0.001). Receiver operating characteristics curve for survival predictive model indicated that the NMELD score is more 
efficient compared with the MELD score in predicting the risk of mortality. Since serum creatinine was not significant in NMELD score, 
further studies to clarify this issue are suggested. 
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1. Background 

Gastrointestinal and liver diseases are among the most 
common causes of morbidity in Iran and constitute a 
substantial proportion of mortality, which imposes enor- 
mous economic consequences (1, 2). Prognosis is an im- 
portant part of the baseline assessment of any disease. 
It is not only the basis of the information that a physi- 
cian provides to the patient, but is also the basis for any 
management method. Proper prioritization of patients 
waiting in queue for a liver transplant (LT) is needed 
and many methods have been developed for this pur- 
pose over several decades (3). In 1964, Child and Turcotte 
proposed a prognostic model for estimation of surgical 
risks in patients with advanced liver disease. Pugh et al. 
proposed a revision of this model in 1973. The modified 



Child-Pugh (CP) prognostic index has been widely used 
for risk stratifying of patients with cirrhosis and to as- 
sess the efficacy of beneficial procedures. At present, the 
CP classification is by far the most extensively applied 
system, as it is easy to use at the bedside (4-8). In 2000, 
the Department of Health and Human Services (DHHS) of 
the USA adopted the 'Final Rule'. According to this rule, 
the primary guidelines for allocation of cadaver livers 
for transplantation should be based on medical urgency. 
Over the years, liver allocation policy has evolved from 
prioritizing liver transplant candidates based on their 
physical location (home, hospital or intensive care unit) 
to medical-based criteria (CP score) consistent with these 
guidelines (5,9). The scoring system of the model for end- 
stage liver disease (MELD) has emerged as an excellent 



Implication for health policy/practice/research/medical education: 

Liver Transplantation is the only treatment option for patients with advanced liver disease. However, there are a few studies on the efficacy of MELD score 
among Iranian patients. If prioritization of liver donation is truly optimized, the mortality rate of patients with advanced liver disease will be consider- 
ably decreased. 

Copyright © 2013, Kowsar Corp.; Licensee Kowsar Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http:// 
creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 



Abolghasemi J et al 



Liver Allocation 



predictor of mortality on the waiting list (WL) (10-13). The 
combination of WL mortality risk and post-transplant 
mortality risk assessed by the MELD score and other fac- 
tors can be used to estimate whether candidates are like- 
ly to derive a survival benefit from a LT. Recently, MELD 
score has replaced CP score in the USA for prioritizing do- 
nor liver allocation (14-27). The introduction of the MELD 
system for transplant allocation in the USA resulted in a 
3.5% reduction in WL mortality while early survival of LT 
recipients remained unchanged despite the selection of 
more ill patients for transplantation (28). In 1991, the first 
LT was carried out in Iran. Currently, more than 10 0 LT are 
being carried out annually in several provinces includ- 
ing Tehran and Shiraz. At present, MELD and CP scores are 
widely used to stratify patients for LT in Iran (29-32). 

2. Objectives 

LT is recognized as the only way to treat patients with 
cirrhosis. The aim of this study was to determine the best 
time for LT and the factors influencing these patients' 
mortality. Optimal allocation of the few available donor 
livers for WL is essential. In this study, the survival rate 
of patients awaiting LT for one, two and three years, and 
the related influential factors were examined. Here, we 
attempt to provide a more efficient scoring system than 
the previous for predicting the mortality risk of patients 
with advanced liver disease. 

3. Patients and Methods 

We evaluated all of the patients with advanced liver dis- 
ease on the WL for LT in Imam Khomeini Hospital (Tehran, 
Iran) during May 2008 to May 2009. They were followed 
up for 3 years. The required data included demographic 
features, cause of cirrhosis and laboratory test results. 
Child-Pugh (CP) score was evaluated by five parameters: 
ascites, encephalopathy, bilirubin, prothomobin time 
and albumin. CP scores of all participants were calcu- 
lated according to the method summarized in Table 1. CP 
scores ranged from 5 to 15 and allowed the categorization 
of the patients into 3 groups (23): a) CP < 7, b) 6 < CP < 9 
and c) CP > 8. 



Table 1. Child-Pugh Scores' Calculating Method 



CP risk class points 


Criteria 


A(l 


B(2 


C(3 




point) 


points) 


points) 


Ascites 


None 


Light 


Large 


Serum bilirubin, mg/dL 


<30 


30-50 


>50 


Serum albumin, g/dL 


>35 


28-35 


<28 


Prothrombin index 


>54 


44-54 


<44 


Encephalopathy 


None 


Minimal 


Advanced 


Total score 


5-6 


7-9 


10-15 



MELD calculations were done on the natural logarithms 
(In) of three variables in blood serum: total bilirubin, cre- 
atinine and the international normalized ratio (INR) of 
prothrombin. MELD score was calculated by the follow- 
ing formula: 

MELD = 9.6 x In (creatinin mg/dL) + 11.2 x In (INR) + 3.8 x 
In (bilirubin mg/dL) + 6.43 

According to the United Network for Organ Sharing 
(UNOS) modifications, the laboratory values below 1.0 
were rounded to 1 to avoid negative scores, and the maxi- 
mum serum creatinine considered within the MELD 
equation was 4.0 mg/dL. In this paper, CP and MELD 
scores were calculated at the beginning of the study for 
all patients. Censorship group included the patients that 
underwent LT during this period or did not die by the 
end of the study. The univariate Cox regression models 
were applied separately to discriminate each significant 
factor. The adjusted survival and hazard functions were 
estimated using a multivariate Cox regression model. 
Factors associated with P-values of less than 0.20 for the 
previous stage (univariate Cox models) were candidates 
for multivariate analysis. In all models, factors with P 
Values < 0.05 were considered to be significant (33). The 
resulting predictive formula was normalized to the same 
scale as the MELD score by linear regression (34). To com- 
pare the new scores with the MELD scores, we used ROC 
curves for survival predictive accuracy. For a binary dis- 
ease outcome, receiver operating characteristics curve 
(ROC) is a popular method for displaying sensitivity and 
specificity of a continuous diagnostic marker. However, 
some disease outcomes are time dependent, and ROC 
curves that vary as a function of time may be more ap- 
propriate. Standard Cox proportional hazard output can 
be used to obtain estimates of time dependent sensitivity 
and specificity, and time dependent ROC curves (35-37). In 
our study, these special Roc curves were used for predic- 
tive accuracy of survival models. About one-third of the 
patients (101), termed group 1, were randomly selected 
for score validation and the remaining patients (204) 
were named as Group 2. Group 2 was used for modeling 
of NMELD scores, and Group 1 was used for its validation. 

4. Results 

Of the 305 patients, 126 (41.3%) were females and 179 
(58.7%) were males. The mean age of the patients was 
40.67 (+ 14.39) (Mean + (Standard Deviation) (M + SD) 
years. The age range was between 18 and 71 years. During 
the study, 71 (23.3%) patients died due to complications of 
liver cirrhosis and 43 (14.1%) underwent LT. Survival rate 
at one, two and three years after enrolling in the WL was 
82%, 73% and 66%, respectively. Survival rate at one, two 
and three years for patients with MELD scores < 10 was 
95%, 90% and 79%, for 9 < MELD < 20 was 89%, 79% and 75%, 
and for MELD > 19 was 55%, 40% and 27%, respectively. 
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Figure 1 indicates that the mortality rate of patients 
with high MELD and CP scores was higher than those 
with lower MELD and CP scores. Using CP scores, 23.9% 
of patients were at stage A, 47.9% at stage B and 28.2% at 
stage C. Wilcoxon test showed no significant difference 
between MELD and CP scores (P = 0.112). In this study, the 
most frequent etiology of cirrhosis was hepatitis, and 
the rate of liver failure associated with viral hepatitis, 
autoimmune and cryptogenic was 34.7%, 16.4% and 22.6%, 
respectively. The detailed data of demographic variables 
are given in Table 2. 



Figure l. Kaplan-Meier Survival Curves for CP and MELD Scores 
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a) Categories of CP scores: l) CP < 6 as class A, 2) CP 7-9 as class B and 3) CP > 
10 as class C b) Categories of ME1D score: l) ME1D < 10, 2) ME1D 10-19 and 3) 
ME1D>20 



Table 2. Demographic, Clinical, Biochemical Features of Pa- 
tients Awaiting Liver Transplantation 



Demographic 



Total 



Group 1 Group 2 



Gender, No. (%) 

Male 180 (58.7) 

Female 125 (41.3) 

Education, No. (%) 

Primary 94(30.8) 



59(58.4) 121(59.3) 
42(41.6) 83(40.7) 

36(35.7) 58(28.4) 



Secondary 

High educa- 
tion 

Marriage status, 
No.(%) 

Single 

Married 

Divorced or 
widow 

Family history of 
the liver disease, 
No.(%) 

Yes 

No 

Cause of cirrhosis, 
No.(%) 

Hepatitis B 
virus 

Hepatitis C 
virus 

Auto Immune 

Cryptogenic 

Alcohol 

PSC 

PBC 

Others 

Clinical Ascites, 
No.(%) 

None 

Light 

Large 

Encephalopathy, 

No.(%) 

None 

Minimal 

Advanced 

Biochemical, 
mean + SD 

Serum biliru- 
bin 

Serum creati- 
nine 

Serum albu- 
min 

INR for prothrom- 
bin time 



168(55.1) 
43 (14.1) 



83(27.2) 
216 (70.8) 
6(2.0) 



52(51.5) 116(56.9) 
13 (12.8) 30 (14.7) 



36 (11.8) 
269(88.2) 



66 (21.6) 

40 (13.1) 

50 (16.4) 
69 (22.6) 
2(0.7) 
27(8.9) 
12 (3.9) 
39 (12.8) 



201 

93 

11 



137(44.9) 
96 (31.5) 
72 (23.6) 



24 (23.8) 
75 (74.2) 
2(2.0) 



19 (18.8) 

17 (16.8) 

1(1.0) 

5(5.0) 

3(3.0) 

15(14.8) 



69(68.3) 
26(25.8) 
6(5.9) 

58(57.4) 
38(37.6) 
5(5.0) 



59 (28.9) 
141(69.1) 
4(2.0) 



17(16.8) 19(9.3) 
84(83.2) 185(90.7) 



27(26.7) 39(19.1) 



14(13.9) 26(12.7) 



31(15.2) 
52 (25.5) 
1(0.5) 
22 (10.8) 
9(4.4) 
24 (11.8) 



132 (64.7) 

67(32.8) 

5(2.5) 

81 (39.7) 
58 (28.4) 
65 (31.9) 



4.290 + 6.394 4.200+ 4.466 + .587 
.480 

0.9 90+1.177 1.024 + 0.924 ±.056 
.101 



3.686 + 2.285 3.617 d 
.123 

1.695 + .744 1.661 d 
.054 



3.821 + .318 
1.761 + .074 



In Group 2, 83 (40.7%) were females and 121 (59.3%) were 
males. 48 (23.5%) patients in this group died due to com- 
plications of liver cirrhosis and 27 (13.2%) underwent LT. 
The mean age of the patients was 41.6 (+ 0.98) (Mean + 
(Standard Error) (M + (SE)) years. The age range was be- 
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tween 18 and 65 years. In univariate analysis, continuous 
variables showed a significant association with the pa- 
tients' age (P = 0.009), In (bilirubin) (P < 0.001), In (albu- 
min) (P < 0.001) and In (INR) (P < 0.001). Categorical vari- 
ables revealed an association between ascites (P < 0.045) 
and encephapolathy (P = 0.013). The results of univariate 
analysis are shown in Table 3. 

Table 3. Univariate Analysis of Risk Factors Associated With 
Mortality in Cirrhotic Patients 



Variables 


Regression 
coefficient 


Regression 
coefficient 
standard 
error 


P 

value 


Ascites (+) 


0.496 


0.242 


0.045 


Encephalopa- 
thy(+) 


1.445 


0.582 


0.013 


Age 


0.014 


0.006 


0.009 


In(bilirubin) 


0.818 


0.131 


< 0.001 


In(albomin) 


-1.836 


0.469 


< 0.001 


In(INR) 


1.555 


0.341 


< 0.001 



Cox multiple regression analysis indicated that there is 
a statistically significant correlation between In (biliru- 
bin), In (albumin) and age, and the risk of mortality. The 
results of multivariate analysis are shown in Table 4. 



Table 4. Multivariate Analysis of Risk Factors Associated With 
Mortality in Cirrhotic Patients 


Variables 


Regres- 
sion coef- 
ficient 


Regression 
coefficient stan- 
dard error 


P- value 


In (bilirubin) 


0.979 


0.153 


< 0.001 


In(albomin) 


-1.745 


0.469 


< 0.001 


Age 


0.033 


0.011 


0.010 



The score of the optimal model was calculated via the 
below formula: 

A = 0.979 x In (bilirubin) - 1.745 x In (albumin) + 0.033 
x age 

By linear regression, we found that the best linear fit 
between A (as independent variable) and the MELD score 
(as dependent variable) was provided by a line with the 
slope of 7.519 and an intercept of 17.661. The correlation 
between the two was good with R= 0.560. Based on the 
above analysis, the scores obtained from the model can 
be derived from the below formula. Here, we've termed it 
the 'New MELD (NMELD)' score: 

NMELD = 10 x (0.736 x In (bilirubin) - 1.312 x In (albumin) 
+ 0.025 x age + 1.776) 

To avoid negative scores, serum bilirubin values below 
1.0 were rounded to 1 and serum albumin values over 5.4 
were rounded to 5.4; the obtained scores were rounded to 
the nearest integer. All analyses for modeling of NMELD 
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Figure 2. Accuracy of the NMELD Score (Dash Line) Using the Co- 
variates of In (Bilirubin), In (Albumin) and age vs. MELD (Solid Line) 
Score Using the Covariates of In (Bilirubin), In (INR) and In (Creati- 
nine). Lines Plot the Estimates of Incident/Dynamic AUC (t) Versus 
Time Under the Assumption of Proportional Hazards. 
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Figure 3. Comparison of the Area Under ROC Curves (AUC) for Pre- 
dicting the Risk of Mortality at3 (AUC = 0.916), 6 (AUC = 0.811), 9 (AUC 
= 0.788) and 12 (AUC = 0.780 ) Months for NMELD Scores. 

were done based on the data collected from Group 2. 
In Group 1, 42 (41.6%) individuals were females and 59 
(58.4%) were males. 23 (22.8%) patients in this group died 
due to complications of liver cirrhosis and 16 (15.8%) un- 
derwent LT. The mean age of the patients was 40.2 (+ 1.53) 
years. The age range was between 18 and 71 years. 

Comparison of accuracy of the NMELD Score (Dash Line) 
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Using the Covariates of In (Bilirubin), In (Albumin) and 
age vs. MELD Score (Solid Line) Using the Covariates of In 
(Bilirubin), In (INR) and In (Creatinine). Lines Plot the Es- 
timates of Incident/Dynamic AUC (t) Versus Time Under 
the Assumption of Proportional Hazards. Figure 2 shows 
that the NMELD score is superior in predicting the risk of 



mortality. Since the AUC values for MELD score are about 
0.5, this model is inefficient for predicting the risk of mor- 
tality. Figure 3 shows that NMELD score is more accurate 
than MELD score in assessing the risk of mortality in the 
short term. Comparisons of the NMELD and MELD scores 
were performed based on data collected from Group 1. 




Figure 4. Comparison of the Area Under the ROC Curves (AUC) for Predicting the Risk of Mortality at 3 (A), 6 (B), 9 (C) and 12 (D) Months Between 
NMELD (Dash Line) and MELD (Solid Line). AUC for Prediction of the Risk of Mortality at 3 (AUC = 0.916), 6 (AUC = 0.811), 9 (AUC = 0.788) and 12 
(AUC = 0.780) Months for the NMELD Score and AUC for Prediction of the Risk of Mortality at 3 (AUC = 0.636), 6 (AUC = 0.697), 9 (AUC = 0.672) 
and 12 (AUC = 0.640) Months for MELD Score. 



5. Discussion 

LT has been accepted as an effective method for the re- 
covery of health in patients with advanced liver disease. 
However, the increasing discrepancy between the num- 
ber of patients on WLs and the number of available donor 
livers has a major impact on the mortality of these pa- 
tients (38). Disease severity at the time of listing has been 
demonstrated as an important predictor of WL death. 



Prognostic assessment of patients with liver cirrhosis is a 
vital subject that often challenges the clinicians. CP score 
is by far the most extensively used both in clinical prac- 
tice and clinical research, and has stood the test of time 
for nearly 40 years. Recently, the MELD score has replaced 
the CP score in the United States for prioritizing liver do- 
nor allocation. Como et al. (39) showed that the use of the 
MELD score produced an advantage for Hepatocellular 
Carcinoma (HCC), but about one in every 11 viral hepatitis 
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patients may be harmed using this scoring system. 

Among the 305 patients in this study, who were listed 
for LT from May 2009 to May 2012, the most common eti- 
ology of liver disease was viral hepatitis (34.7%). This is in 
contrast with the findings of the Malinchoc et al. study, 
conducted at a medical center in the USA where alcoholic 
liver cirrhosis was the most prevalent (38.8% to 90.5%) fol- 
lowed by viral hepatitis (4.8% to 10.4%), and also with the 
results of a study by Sumskiene et al. in Lithuania where 
alcoholic liver cirrhosis was the most prevalent (28.9%) 
followed by viral hepatitis (17.8%) (40, 41). This major dif- 
ference is mainly due to the prohibition of alcohol by 
Islamic rules; therefore, alcoholic liver cirrhosis is not 
common in Iran. According to the data shown in Figure 4, 
NMELD score is superior to the MELD score in accurately 
predicting the risk of mortality in Iranian patients with 
advanced liver disease. 
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